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We calculate the spin relaxation time of conduction electrons in n-doped bulk gal- 
lium arsenide. We consider the Elliot-Yafet spin-relaxation mechanism, driven by 
Coulombic-impurity and electron- electron scattering. We find that these two scat- 
tering mechanisms result in relaxation times of equal order of magnitude, but with 
disimilar dependences on doping density and temperature. Our theoretical results are 
compared with experimentally measured spin relaxation times in gallium arsenide. 

PACS numbers: 72.25.Rb, 76.30.Pk, 71. 10. Ay, 03. 67. -a 

There is currently a great interest in the properties of semiconductors derived 
from the electronic spin. This interest stems largely from the foreseen growth of 
the scope of applications of the spin degree of freedom in electronics (spintronics) 
and computer science (quantum computing) p!J |2j. Among the properties of in- 
terest in relation to potential applications, a central place is occupied by the spin 
relaxation time, i.e., the characteristic time a spin-density imbalance lasts inside a 
given material or structure. Clearly, long relaxation times are generally desirable, 
and an important progress in the search of systems with long relaxation times has 
been made in a series of recent experiments [3J IU E|. In those experiments it was 
found, among other things, that the spin-relaxation time in semiconductors can be 
extended by more than two orders of magnitude with appropriate negative-type 
doping. 
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Several authors have recently analyzed this issue from a theoretical point of view. 
Song and Kim [5| reviewed the known spin-relaxation mechanisms and constructed 
a phase diagram to depict graphically which mechanism is dominant in each portion 
of parameter space. According to those authors, the D'yakonov-Perel (DP) [7j mech- 
anism dominates over the Elliot- Yafet (EY) one [HI Ej (driven by electron-impurity 
scattering) except at low temperature in n-doped zincblende semiconductors. Wu 
and Ning (TO) studied D'yakonov-Perel spin relaxation in n-doped GaAs. Their main 
finding, as far as explaining the above-mentioned experiments is concerned, is that 
the DP mechanism produces and increasing relaxation time as a function of applied 
magnetic field, in qualitative agreement with experiment at high density [3]. A re- 
lated study, based on a kinetic equation formalism and including the k-dependence 
on the conduction-band g-factor, was done by Bronold et al. [TT]. These authors find 
a good agreement between their theory and experiment [3] for the relaxation time as 
a function of magnetic field at high doping density. While the high-doping-density 
case can be treated in terms of delocalized electrons — electrons at the bottom of the 
conduction band — , the low-density regime calls for a description based on electrons 
localized at the impurity sites. Kavokin and Gor'kov and Krotkov ^Hj studied 
this regime. The latter authors conclude that the Dzyaloshinskii-Moriya interac- 
tion between localized electrons does not explain the experimentally observed spin 
relaxation times, as was claimed by Kavokin. 

Here we report our calculations of the spin relaxation time of electrons in the 
conduction band of n-type semiconductors caused by the Elliot- Yafet (EY) scat- 
tering mechanism [HI Ej mediated by both electron-impurity and electron-electron 
interactions. Our calculations extend those of Chazalviel [Hj and Boguslawski JT3] 
for electron-impurity and electron-electron scattering, respectively. Both of these 
authors considered only the case of almost equal up-spin and down-spin popula- 
tions. Furthermore, Chazalviel studied only the zero temperature limit while Bo- 
guslawski treated the high temperature case by introducing Boltzmann rather than 
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Fermi distributions. We generalize those calculations by considering arbitrary spin 
populations described by finite temperature Fermi distributions. Also, while those 
authors were mainly interested in indium antimonide, we concentrate here on gal- 
lium arsenide, which is one of the preferred semiconductor alloys in recent studies 
of spin relaxation. 

As it is clear that only an incomplete understanding of the experiments reported 
in Ref. [I] has emerged so far, our aim here is to provide a missing piece in the 
set of theoretical scenarios considered up to now. We point out that none of the 
recent theoretical studies of spin relaxation in n-type GaAs has taken into account 
the electron- electron EY mechanism, which had been studied by Boguslawski ^S] i n 
a different context. 

The Coulomb interaction is independent of spin and therefore cannot cause spin- 
flip transitions between conduction-band states that are spin eigenstates. The pos- 
sibility of a spin flip in the EY scattering mechanism arises from the fact that the 
conduction-band states of some semiconductors are not spin eigenstates, which in 
turn is due to the spin-orbit contribution to the crystal Hamiltonian. The conduc- 
tion band states of zincblende semiconductors are linear combinations of both spin 
eigenstates with coefficients that are functions of the crystal momentum k. Usu- 
ally one of the two components in these admixtures has a much larger amplitude 
than the other. Therefore the mixed states retain a clear correspondence to the 
original pure-spin states, and are still referred to as the "spin-up" and "spin-down" 
states. We denote these states by |k±). The matrix element of a spin-independent 
scattering potential (like that of ionized impurities and electron-electron interac- 
tion) between these states can thus be non-zero, thereby activating the spin-flip EY 
scattering mechanism. 

Our general goal is to understand the spin relaxation process in n-doped bulk 
semiconductors studied recently in ultrafast Faraday rotation experiments [4]. In 
this article we do not investigate the initial stages of the relaxation process, when 
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electron- hole scattering is predominant, but concentrate rather on the slower decay 
that takes place after the electron-hole recombination has been completed. There- 
fore, the starting point of our calculation is the assumption that there is a doped 
unpolarized electron gas of density n present in the host semiconductor, and an 
additional density n ex of photoexcited electrons occupying conduction-band states 
of the type |k+). Since we "look" at the system at least several hundreds of pi- 
coseconds after the pump pulse has created the photoexcited polarized electrons, 
we assume that the two spin species have reached thermal equilibrium and can be 
described by Fermi-Dirac distributions /±(k) with a common temperature T. We 
restrict our treatment to cases where n ex < n. 

We are interested in the relaxation of a spin-density difference = n + — n_, 
where n± = (1/V) Z)k/±(^) are the densities of spin-up and spin-down electrons 
in the conduction band. We assume that initially = n ex . The relaxation rate is 
defined by 

1 _ fid h+ — h- 2n + 
Ti n d n d n d ' 

Here we have used that n + + n_=0 since no further excitation nor recombination 

take place in the considered time regime. 

For electron-impurity scattering, the relaxation rate is calculated as 

1 2n ■ 

t = — -EE/-(^)[i-/ + W]W-/ + W[i-/-(fc')]W 

H ri d k k , 

= ^EEl/+(*)-/-(*0]w, (2) 

71(1 k k' 

where we have denoted the density of impurity scatterers by rij, and used that 

Tk' ,k+ — ?k+- +k'-- The spin-flip transition rate due to the impurity EY mechanism 

is 

2-7T 9tt 

T k+ ^k'- = i- 5(E k ,-E k ) |(k / -|y i |k+)| 2 « — 5(Ev-E k ) ^(k-k') 2 |(k'-|k+)| 2 
n Ti 

(3) 

The electron-impurity scattering is caused by the Coulomb potential of the ionized 
silicon donors 16J , whose Fourier transform is given by Vi(q) = ATre 2 /eV (q 2 + k 2 ), 
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where e is the lattice dielectric constant, V is the volume, and the screening wave 
vector in the quantum Debye-Hiickel theory for a degenerate electron gas is k s = 

1 /2 

(67r e 2 n/e Ep) ■ For n and Ep in this expression we take the values corresponding 
to the doped unpolarized electron gas. The spin-mixed conduction-band states of 
zincblende semiconductors can be calculated with the kp perturbation theory [17J. 
The scalar product of opposite-spin states is given by [TU [TH] 

(V-\k + ) = ^-(k z k> + -k> z k + ), (4) 

where 7 = 2A(A + 2£ G )/(A + E G )(2A + 3E G ), k + = k x + ik y , m* is the conduction- 
band effective mass, A is the valence-band spin-orbit splitting, and E g is the bandgap 
energy. Substituting Eq. (jlj into Eq. (J3J) yields 



T k+ ^- = ^5(^-E k )^(k-k') 2 (k 2 k' z 2 +k' 2 k 2 z -2k z k' z k-k'). (5) 



Doing the integral over k' and the angular part of the integral over k one obtains 



1 rti hj 2 e 4 



dkk'\f + (k) - f-(k)] 



-2 + (1 + a) In 



2 + a 



(6) 



T\ rid 67rm*E G e 2 

where a = (k s /k) 2 /2. 

For electron-electron scattering, the relaxation rate is calculated from 

^ £{/_(k)/ CT (k') [1 - / + (k + q)] [1 - / CT (k' - q)] T k _ ; k v-, k+q+; 

- / + (k)/ ff (kO [1 - /_(k + q)] [1 - / CT (k' - q)] T k+; k ^k +q - ; k'- q J , (7) 



dt 

kk' 



where T kl ± . k ' i(J ^ k2 T ; y a- are the transition probabilities of scattering events where 
one of the electrons flips spin and the other does not. The likelyhood of double 
spin flips can be neglected. Taking into account the direct and exchange terms, the 
transition probability according to Fermi's golden rule is given by 

27T 

T\t- ■ k',cr^k+ q ,+ ;k'- q ,cr = -^-&{E f — Ei) ■ 

■ |(k+q, +; k'-q, a\ V |k, k', a) - (k'-q, a; k+q, +| V|k, -; k', a)\ 2 . (8) 
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The interaction potential here is again the screened Coulomb potential given above. 
Its matrix elements between two-particle states are given again by Eq. (j3J) after a 
few standard approximations are made ^H]. The electron-electron EY relaxation 
time is obtained by solving Eq. (JJJ) numerically [TH] . 

We first present results for the impurity EY mechanism obtained numerically 
from Eq. (JBJ. That equation is valid for all zincblende semiconductor alloys, but in 
this article we concentrate on gallium arsenide, given the relative importance of this 
material in recent experimental studies. For a given material, the impurity EY spin 
relaxation time depends on the dopant density Hi, the density of the excess electrons 
with spin up, = n ex , and the temperature T. The density of donated electrons, 
n, which also enters implicitly in the calculation, is taken to be n = rij. We envision 
that the initial spin-density difference, rid, comes from the photo-excitation caused 
by a circularly-polarized pump laser pulse in a pump-and-probe experiment. 

Our calculations show that, for impurity EY scattering, the dominant factors in 
the determination of T\ (according to (|H)l) are the densities n and ?v Therefore, we 
show first the time T\ versus n, for three different values of rid and of T (Fig. (1)). 
At low temperature and low spin-density imbalance, T\ decreases by a factor of 1000 
for a variation of a factor of 100 in n; thus, T\ follows roughly a power law on n 
with an exponent of -3/2. This dependence comes from both n and rij. The density 
of scattering centers, rii, enters linearly in Eq. (JB|L while the density of donated 
electrons n affects the integral in that equation through the screening length l/k s 
and the Fermi functions. 

The dependence of T\ on T is strongest at low values of n, as could be expected 
from reasoning on the effect of temperature on Fermi functions for different densities. 
For n = 10 18 cm -3 there is practically no temperature dependence up to the highest 
temperature considered here, T =100 K. The relaxation time also depends very 
weakly on rid, except when rid approaches n. 

In Fig. 2 we present more detailed data on T\ versus the spin-density imbalance, 
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rirf, for doping density n = 10 16 cm~ 3 . Temperatures from 1 K to 100 K are con- 
sidered, as indicated in the figure caption. The main conclusion we can draw from 
Fig. 2 is that there is essentially no variation of T% with rid for rid < n/10. (This 
conclusion is also valid for n = 10 17 , 10 18 cm -3 , data not shown.) The variation of 
n between rid — n/10 and rid — n is roughly of a factor 2 at low temperature, and 
almost none at high temperature. (The definition of low and high T depends on n, 
as mentioned above.) 

We shall now move on to the results for the electron-electron EY scattering. 
We define the relative density n exr as the ratio between the density of photoexcited 
electrons and the doping density, n exr = n ex /n. Figure 3 shows the spin relaxation 
time Ti, obtained by solving Eq. (J7J) and using Eq. ((TJ), as a function of n for 
n exT = 0.001,0.01,0.1 and various temperatures between 1 K and 315 K. Once 
again we obtain that T\ decreases rapidly with doping density much as it does in 
the case of impurity EY scattering. On the other hand, Fig. 3 shows also a strong 
temperature dependence, not present for impurity scattering. 

In Fig. 3 we can identify different regimes that can be broadly classified according 
to the doping density and the temperature. A high temperature regime appears for 
T > 10 K and n<10 16 cm -3 , or T > Tp, where Tp is the Fermi energy of the 
doped electron gas. The low temperature regime is located on the "opposite corner 
(upper-right)" of the plot, and is given roughly by T < 3.2K and n>10 17 cm -3 
(T < Tp). Other combinations of n and T correspond to T w Tp. The following 
main conclusions can be drawn from Fig. 3: (i) At high temperature (Tp < T), 
T\ does not depend significantly on n exr ; (ii) At low temperature (T < Tp), as n 
increases T\ decreases and this decrease is sensibly greater the higher the value of 
n exr . (iii) At low temperature (see curves for T = 1 K, 3.16 K), high 

Tlexr (see 

riexr = 0.1), and high density (n ~ 10 18 cm~ 3 , the temperature dependence of T\ 
dissapears. We have analyzed in detail the dependence of T\ on n, T, and n exr , but 
at this point we shall just refer the interested reader to Ref. |18j . 
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Finally, we compare our theoretical results with the corresponding experimental 
findings of Ref. The longest spin- relaxation times reported in Ref. 0] where 
obtained at low magnetic field and with a doping density of n = 10 16 cm -3 . In Fig. 4 
we plot the relaxation time versus temperature. The two solid curves are results from 
our electron-electron EY scattering calculation and the triangles are experimental 
data from Ref. jlj. There is a large disagreement between the experimental points 
and the calculated curves for all temperatures, of two to three orders of magnitude. 
We recall that the impurity-EY relaxation times (Figs. 1 and 2) are of the same 
order of magnitude as the electron-electron EY relaxation times shown in Fig. 4, 
and therefore they also compare unfavorably with experiment. The fact that both 
types of EY relaxation times are much longer than the ones measured experimentally 
could indicate that another, more efficient, spin-flip mechanism dominates the spin 
relaxation in this system. A partial agreement between experiment and theory has 
been found in terms of the D'yakonov-Perel spin-flip mechanism in a calculation that 
includes a careful treatment of the semiconductor band structure and experimentally 
measured mobilities as input parameters ^H]. That treatment seems to fit the 
experimental relaxation times at a temperature of around 100K, but is off by two 
orders of magnitude at low temperature. Therefore, we cannot rule out yet that a 
refined version of the EY calculation presented here could help explain the measured 
relaxation times in n-doped GaAs, particularly at low temperature. The refinement 
in the calculation should probably be made on the expression of the spin-mixed 
conduction-band states, the treatment of screening of the Coulomb potential, and 
also the Born approximation to scattering rates used here and in Refs. [T4^ ITH] . 
Furthermore, at low density the possibility of electron localization at impurity sites 
has to be taken into account ^21 EH! • 

In Ref. f 1 it is mentioned that the electron-electron EY mechanism could help 
understanding the temperature dependence of the relaxation time at low temper- 
ature (below 30 K). However, as noted by the authors, the estimates based on 
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the treatment of Boguslawski which assumes classical statistics, yield a re- 
laxation rate that is actually too fast compared to experiment. Our calculation 
of the electron-electron EY mechanism shows that Boguslawski's results cannot be 
used to interpret experiments at low temperatures (< 100K). Indeed, the use of 
Fermi distributions enforces Pauli blocking which greatly reduces the efficiency of 
electron-electron scattering as compared to the case of Boltzmann statistics. 

In conclusion, we have generalized the treatments of Chazalviel ^3] and Bo- 
guslawski of impurity and electron-electron EY spin-flip scattering in n-doped 
zincblende semiconductors. Our treatment of spin relaxation takes into account the 
dependence of T\ on doping density, spin-density imbalance, and temperature. We 
apply this theory to GaAs and analyze the main features of the dependence of T\ 
on those system variables. This theory should mainly be considered as a starting 
point for more refined treatments of the Elliot- Yafet mechanism, as a complete un- 
derstanding of recent experiments [3J E] based on the known spin-flip mechanisms 
in semiconductors has not emerged so far. 
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FIGURE CAPTIONS 



FIGURE 1. Impurity EY spin relaxation time of GaAs versus doping density at three 
different temperatures (5 K, 50 K, and 100 K) and three spin-density differences: 
rid — 10 14 cm -3 (solid lines), 10 15 cm -3 (dotted lines), and 10 16 cm -3 (dashed lines). 

FIGURE 2. Impurity EY spin relaxation rate of GaAs versus spin-density difference 
{rid = Uex) for doping density n = 10 16 cm -3 (equal to the impurity density rii.) The 
various curves correspond, from top to bottom, to T = 1, 10, 20, ... , 100 K. 

FIGURE 3. Electron-electron EY relaxation time versus doping density n for several 
values of the temperature and the relative density of photoexcited electrons, n exr . 

FIGURE 4. Spin relaxation time versus temperature. Comparison between experi- 
mental data (Ref. jl]) and results of the electron-electron EY scattering calculation. 
The theoretical curves are for n = 10 16 cm -3 and n er = 10 14 and 10 15 cm -3 . 
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